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I present recent progress in theoretical modelling of cosmological density-velocity relations in 
£ — ■ the weakly nonlinear regime. The relations are local, based on rigorous perturbation theory 

and include the effects of smoothing of the density and the velocity fields. For small smoothing 
scales, they can be improved by slight adjustments based on N-body results. The relations 
can be useful for density-velocity comparisons, yielding a fair estimate of £1 and offering a 
method for disentangling Q and bias. 



1 Introduction 
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In the gravitational instability paradigm for the formation of structure in the Universe, the pecu- 
liar motions of galaxies are tightly related to the large-scale mass distribution. The comparison 
between the density and the velocity fields can serve as a test of the gravitational instability 
hypothesis and as a method for estimating the cosmological parameter Q (Dekel et al. 0) . In 
linear theory, the relation between the density and the velocity fields is 

I" ^(x) = -/(Q)- 1 V-v(x), (1) 

where /(O) ~ fi 0-6 . The relation between the density contrast and the velocity divergence is 
thus linear and local. 

This equation is applicable only when the density fluctuations are small compared to unity, 
°S = (5 2 ) <C 1. However, sampling of galaxies in current redshift surveys and random errors 
in the peculiar velocity data enable reliable dynamical analysis at smoothing scales down to a 
few Mpc, where fluctuations already exceed the regime of applicability of linear theory. On 
the other hand, even at scales as small as a few /i" 1 Mpc, 'typical' (rms) fluctuations are still 
not in excess of unity|] One could thus hope that to construct the relation between the fields in 
question, perturbation theory can be effectively used. 

In perturbation theory, one approximates the evolved density contrast as a sum of terms Sj , 
each corresponding to the j th order in perturbation theory, 

6 = 5i + 6 2 + 5 3 + ■■■ ■ (2) 
Similarly, defining a variable 9 proportional to the velocity divergence, 

e = -/(n)- 1 v-v, (3) 

we can expand it in a perturbative series, 

9 = 9 1 + 9 2 + 9 3 + . . . . (4) 

"For a Gaussian smoothing function, commonly used in the density-velocity comparisons. 



The solutions for Sj and 6j were obtained more than a decade ago (Fryi, Goroff et al. i). In 
general, the j order solution is of the order of ~ a\. In the weakly nonlinear regime, i.e., 
when as < 1, which is our case, we may hope the perturbative series to converge rapidly. 

Linear theory is nothing but perturbation theory truncated at the lowest, i.e. the first order. 
A natural way of extending the linear relation (|l]) is to include higher-order contributions to 
the density contrast and the velocity divergence. Unfortunately, these contributions are non- 
local, hence the nonlinear density vs velocity-divergence relation (DVDR) at a given point has 
a scatter. Therefore, as a local estimator of the density from the velocity divergence we adopt 
the conditional mean, i.e., the mean density given the velocity divergence. 

The calculation of the conditional moments is a statistical problem. To compute mean 5 given 
6, (S)\g, we need the joint probability distribution function (PDF) p(5, 6). Under an assumption 
of Gaussian initial conditions, adopted in this work, this function is initially a bivariate Gaussian 
(degenerated, since Si = 9\). Subsequent nonlinear gravitational evolution, however, drives the 
PDF away from its initial shape. Our problem here is to calculate it in the weakly nonlinear 
regime. 

Mathematically, a PDF is given by the inverse Fourier transform of its characteristic function 
which in turn is related to the cumulant generating function K, by the equation <I> = exp [K]. 
The cumulants, K mn , from which IC is constructed, 
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E ^r^TiisT, (5) 

lll.lt. 

(m,n)^(0,0) 

are the connected (reduced) part of the joint moments of S and 0, K mn = {S m n ^ conn . We 
introduce K, because in the weakly nonlinear regime, the cumulants obey a scaling hierarchy in 
as such that the series (|5|) is a power series in as (even in the case of the standard variables 
\i = 5/ as and v = 0/gq). Therefore, truncating @ at order p we neglect contributions which 
are 0{a^ +l ). This yields the truncated PDF and in turn the weakly nonlinear DVDR up to a 
given order in perturbation theory. 



2 Results 

Up to third order (Chodorowski & Lokasi) we have 

(S)\ e = a l e + a 2 (e 2 -a 2 e ) + a 3 e 3 , (6) 

where a 2 is the variance of the velocity divergence field. The reverse relation (Chodorowski et 
al.u) is 

(e)\ 5 = r 1 S + r 2 (6 2 -a 2 5 )+r 3 5 3 (7) 

(an extension of a second-order formula by Bernardeau 0) . Because the DVDR has a scatter, 
expression (0) cannot in general be obtained by direct inversion of expression (g). In density- 
density comparisons, one reconstructs the mass density field from the observed velocity field 
and compares it to the observed galaxy density field. In velocity-velocity comparisons, one 
reconstructs (under some assumptions about bias) the velocity field from the observed galaxy 
density field and compares to the observed velocity field. Formula @ is thus relevant for 
density-density comparisons, while formula ([?]) is relevant for velocity-velocity comparisons. 
Indeed, for irrotational flows the velocity field can readily be recovered from its divergence given 
some boundary conditions at large distances, 




Figure 1: Left: joint probability distribution function-ifor the density contrast (8) and the velocity divergence 
(9) from an N-body simulation by Chodorowski et al.,B shown against 1 + 8 and 1.5 + 8. Right: the coefficients 
dj (see eq. ||), as functions of the smoothing scale, from N-body simulations by Chodorowski & Stompor. u 
Simulation results are shown as thin curves; third order perturbation theory predictions are shown as thick lines. 
The simulations are for power-law power spectra with the spectral index n equal to —1.5 (left) and —1 (right). A 

top-hat smoothing function is used. 



The coefficients aj and rj are the previrialization-, skewness-, and kurtosis-like combinations 
of the joint moments of 5 and 9. Calculations of the reduced moments have already become a 
small industry in perturbation theory. For Gaussian smoothing, the related moments for 5 and 
for 9 (separately) were computed by Lokas et a/.0 and Lokas et aZ.0 The coefficients aj and rj, 
though constructed from the joint moments, possess similar mathematical structure and were 
computed in an analogous way. Hi 

Chodorowski et al. a also calculated the scatter in the DVDR and found it to be relatively 
small.]] Nevertheless, one can try to further reduce the variance of density estimators by including 
off-diagonal components of the velocity deformation tensor. Ind eed , including the shear of the 
velocity field reduces the scatter in the density-velocity relation. I 



3 Tests against N-body simulations 

Perturbative expressions for density-velocity relations should be checked against N-body simu- 
lations, in order to assess their range of applicability. The results of the simulations by Chodor- 
owski et al.u and Chodorowski & Stompor □ (see Figure ||) can be summarized as follows: 

• For large smoothing scales (as <C 1), the numerical coefficients aj and r,- converge to the 
predicted values. (In this sense, the predicted DVDRs are 'asymptotically unbiased'.) 

• For smoothing scales as small as a few megaparsecs (as ^ 1), 

— The numerical coefficients slightly, but systematically, deviate from predicted values. 
(Apparently, higher-than third-order contributions become significant.) 

— On the other hand, third-order polynomials, with fitted coefficients, provide excellent 
fits to N-body data. (Higher-order contributions modify the values of the coefficients 
for j = 1,2,3, but do not significantly induce non-zero values for j > 3.) 



The ratio of the scatter to the rms fluctuation of the density field, as, vanishes in the limit as — > 0, as 
expected. 



4 Implications for the value of Q 



For the sake of simplicity, here we will consider only second-order corrections to the linear 
relation (|l|). Neglecting third-order terms in equation (^) yields 

8 = 5 - a 2 (5 2 - a 2 s ) (9) 

(r2 = — ci2, 0.2 < a2 < 0.3 iH). Assuming linear bias between the galaxy density field, S g , and 
the mass density field, 5 g = b5, we can express the relation in terms of the observable quantities 
V • v and 5 g , 

- V • v = l3S g - a 2 b- 1 f35 2 g . (10) 

Here, b) = f{Q)/b and we also neglect the offset in the relation. We see that the quadratic 
correction has the opposite sign to the linear term. It implies that the linear relation overesti- 
mates the predicted velocities. Hence, in order to match the observed ones, a lower value of Bis 
needed. In effect, applying linear theory biases the estimate of Q. low (cf. Strauss k. WillickO). 

Note that the two coefficients of the binomial in 5 g on the RHS of equation (|To| ) are different 
combinations of (3 and b, or /(fi) and b. This offers a way to solve for and b separately! 

5 Summary 

The weakly nonlinear estimators of density from velocity, and vice versa, presented here, 

• are asymptotically unbiased; 

• for smaller smoothing scales, can be improved by slight adjustments from N-body data; 

• should be used in density-density and velocity-velocity comparisons, in order not to bias 
the estimate of low; 

• offer a method for disentangling 0, and bias. 
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